Quantum state stability against decoherence 
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We study the stability of the coherence of a state of a quantum system under the effect of an 
interaction with another quantum system at short time. We find an expression for evaluating the 
order of magnitude of the time scale for the onset of instability as a function of the initial state of 
both involved systems and of the sort of interaction between them. As an application we study the 
spin-boson interaction in the dispersive interaction regime, driven by a classical field. We find, for 
this model, that the behavior of the time scale for the onset of instability, with respect to the boson 
bath temperature, changes depending on the intensity of the classical field. 

PACS numbers: 03.67.-a, 03.65.-w 



I. INTRODUCTION 

The question as to whether or not a pure quantum 
state can persist in the macroscopic world has been con- 
sidered since Schrodinger introduced his gedanken exper- 
iment known as Schrodinger 's cat p| and it was in this 
way that the entanglement concept was introduced into 
the quantum world. Decoherence is a word to indicate 
that the state of a quantum system is not pure. In gen- 
eral, a hamiltonian interaction which generates entangle- 
ment between two quantum systems produces reversible 
decoherence in each of the involved systems. However, 
when a system interacts with an infinite number of sys- 
tems or a single system interacting with an infinite num- 
ber of degrees of freedom, the concept of irreversible deco- 
herence arises, which is usually called decoherence only. 
Thus, decoherence has become the terminology for the 
irreversible evolution of a quantum state due to its inter- 
action with an environment Q- In this context numer- 
ous works allow us to gain intuition about the dynamical 
behavior of open quantum systems. For instance, spon- 
taneous emission arises from the coupling of a system to 
a noise vacuum environment [j| , Dalvit et al. 4] studied 
various measures of classicality of the pointer states of 
open quantum systems subjected to decoherence. Also 
J. I. Cirac et al. [j| found a dark state of a single two- 
level ion trapped in a harmonic potential by controlling 
its motion degree of freedom. This is achieved by gen- 
erating a squeezed motion state. T. A. Costi and R. H. 
McKenzie gave a quantitative description of the en- 
tanglement between a two-level system and an environ- 
ment for an ohmic coupling. D. DiVincenzo and D. Loss 
provided an exact analysis of the weak coupling limit 
of the spin-boson model for an ohmic heat bath in the 
low temperature limit, using non-Mar kovian and Born 
approximations. S. Bose et al. Q studied the enforcer of 
entanglement between a two-level system and a quantized 
mode in a thermal state. Decoherence-free subspaces in 
cavity QED have been found @. Also a method has 
been developed which achieves the slowing down of de- 



coherence and relaxation by fast frequency modulation 
of the system- heat-bath coup ling [lCj. In addition A. G. 
Kofman and G. Kurizki developed a unified theory 
of dynamically suppressed decay and decoherence by an 
external field in qubits to arbitrary thermal bath and de- 
phasing sources. The problem of stability of a quantum 
state under a class of Lindblad dissipative dynamics has 
been also studied [f3 ]. 

In this article we study, in a simple form, the problem 
of coherence stability, at short time, of a quantum state 
without Markoff nor Born approximations. We have 
found an expression which allows us to estimate the time 
scale order of onsetting the decoherence for a quantum 
system when it interacts only with another single quan- 
tum system. The expression can also be applied when 
the second interacting system is in a thermal equilibrium 
or is in an adiabatic dynamics regimen with respect to 
the studied one. We apply our result to the following 
cases: (i) the reversible pure-dephasing interaction with 
a numerical simulation, (n) a cavity mode driven by a 
thermal light, and (Hi) the spin-boson dispersive inter- 
action driven by a resonant classical field. 



II. DECOHERENCE RATE 

The stability against decoherence is understood to be 
the process where quantum coherence is preserved along 
evolution |12|. So, we can say that an initial pure state 
\ip) is stable against decoherence during the time td if 
trp 2 (t) « 1 for all t < td- The decoherence of a state, rep- 
resented by p, is measured by means of the first-order en- 
tropy: s(t) = l—trp 2 (t). In order to consider the stability 
of the coherence, we assume that s(t) is an analytic func- 
tion of t |l2j • To find out the td stability time scale order 
for which s(t) remains being approximately s(0) = 0, it 
can be expanded as a Taylor series, in such a way that 
the time scale order will be given by td — 1/ \/s^(0), 
where s n (0) denotes the nth derivative of s(t) evaluated 
at t = 0, and s„(0) is the lowest order derivative differ- 
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ent from zero. The first two derivatives of the first-order 
entropy are: si = —2tr(pp) and s 2 = —2tr(p 2 + pp), with 
p and p denoting, respectively, the first and the second 
derivatives of p with respect to t at time t. 

Now let us suppose that a system under study, labeled 
by a, is interacting with the R system through the V(t) 
hamiltonian. We consider the whole a — R system t be 
isolated or: R being in a thermal equilibrium or in an 
adiabatic dynamics regimen. Initially the a system is in 
a pure state \ip) and the R system is in the state p R . In 
the whole tensorial product Hilbert space, H a <g> TLr, the 
dynamics of the composite system state p(t) is driven by 
the Master equation (H= 1): 

p(t) = -i[V(t),p(t)}, (1) 

whereas the dynamics of the partial density operator of 
the a system, p a (i) = trn(p(i)), is governed by 

p a (t)=-itr R [V(t),p(t)], (2) 

where tr R denotes the tracing up over the R system. 

To find out if the initial pure state \tp) is stable un- 
der the dynamics described by the and J2J Eqs., and 
besides, to obtain the expression for the td, we consider 
the change in s(t) by calculating its derivatives at t = 0. 
Since we suppose initially the state of the system a to 
be pure, then s(0) = and si(0) = 0; hence the sec- 
ond derivative of the first-order entropy, which could be 
different from zero, is given by 

S2 (o) = 2>{{v% mn) R - (v) R ([v, mm R ) a , (3) 

where we have denoted the average on the system a (R) 
by the subindex a (R), and they are taken at t = 0, 
and V = V(0). Thus, when the \ip) pure state does 
not commute with V(0), e.i. [V, |V>) (V'l] 7^ 0, the value 
of td = 1/ \Zs 2 (0) gives a time scale order for the onset 
of decoherence. We can also see from Eq. that, in 
principle, td © will be a functional of both the initial 
pure state of a as the initial state of the R system, and 
is proportional to 1/g. It is important to point out that 
S2(0), Eq. ©, is the positive defined operator-correlation 
in the R system between the V and [V, IV'XV'I] operators, 
averaged on When the state commutes with 

V(0), one must calculate the third derivative of the first- 
order entropy at t = in order to have the time scale 
order. 

Here it is worth emphasizing that td was found using 
neither Markovian nor Born approximations. It is well to 
recall that the Markov approximation is a coarse grained 
dynamics description, in the sense that the time scale 
in which the system is observed is much longer that the 
characteristic correlation time of the reservoir. Since, 
in that case, fine temporal structure can not be seen, 
si(0) ^ in general. Our description of the dynamics of 
the system corresponds to a time scale smaller than the 
correlation time of the reservoir. 

By way of examples first let us consider familiar in- 
teraction models which generate reversible decoherence, 



that is, a far from the resonance interaction which is 
known as simplest pure-dephasing mechanisms Il5l 
EH described by the V = g(b + &t) a z hamiltonians 1151 . 
where a z is the z-component of the a spin-1/2 operator 
with eigenstates |0) and |1). b and b> are the boson anni- 
hilation and creation operators respectively, and g gives 
account of the effective coupling strength. For this case 
the decoherence time scale order is given by 



2g^((A(b + b^) R ^l-(^\a z \^ 

where yj ((A(6 + fot)) 2 ) r is the root-mean-square devia- 
tion of the b + b' boson quadrature at t = 0. 

We can see that the initial coherence of the boson state 
plays an important role, i.e., for a fixed state, an 
initial boson squeezed state [H, [T]j causes a decoherence 
time scale smaller than one caused by an initial Fock state 
or by a thermal state with equal average boson number. 
In this effective model, |0) and |1) states are affected 
only by a phase and each one is stable under this pure- 
dephasing mechanism. So, from we can see that, for a 
fixed boson state, the states on the equator of the Bloch 
sphere have a smaller decoherence time scale than the one 
for states being near to the poles. It is worth noting that, 
for this pure-dephasing mechanism, any two states, \tp) 
and its orthogonal \ip±), have the same td H(j. For the 
Jaynes-Cummings resonance |T3| interaction model, one 
can show that two states, and its orthogonal \ip±), 
have different td poj ]. 

Fig. shows the exact evolution of the s(t) first- 
order entropy of the two-level system under an effective 
pure-dephasing interaction with a boson field mode. The 
two- level system is initially in the |+) eigenstate of a x 
and the boson mode is in: a Fock state (solid), ther- 
mal state (dash), and vacuum squeezed state (dot), each 
one with the same average boson number (n) = 3. The 
r squeeze parameter is such that yj ((A(b + b^)) 2 ) R = 

ss 0.26795 and for thermal and Fock state 
V /((A(6 + fot))2) jR = ^7 ~ 2.64575. Thus, in this partic- 
ular case the td time scale orders differ in one order of 
magnitude as can be seen in Fig. ^ 

As a second example, let us consider the explicit model 
which consists of a cavity mode driven by a thermal light 
|l2l H^ l . In the Markov limit and Born approximation 
this model accounts for cavity losses. The hamiltonian 
in the interaction picture of this physical model is given 
by 

V(t) = 5^ ft- (art e*^-")* + aV.-e"^ (5) 

3 

where w is the mode frequency and Wj is the frequency 
of the jth mode of the thermal light, a and a* are the 
annihilation and creation operators respectively, and rj 
and rj are the annihilation and creation operators respec- 
tively of the jth mode of the thermal light, gj account for 
the effective coupling strength between the main mode 
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FIG. 1: Evolution of s(t) for the two-level system under an 
effective pure-dephasing interaction with a boson field mode. 
The initial state of the boson mode is: a Fock state (solid), 
a thermal state (dash), and a vacuum squeezed state with 
quantum noise reduced in the b + 6* quadrature. 

and the jth mode of the thermal light. For this model 
the decoherence time scale is 



2 v /(7 + 2 7 T)«a t a) - (at)( a ))+ 7T ' 

being 7 = ^ . \g } -| 2 which gives account of the whole mag- 
nitude of effective couplings strength, jt = J2j l3j| 2 ( n .i ) > 
and they are the reservoir correlation functions |13j . 
From Eq. © we can see that, at zero temperature 
((rij) — Vj), a coherent state is stable at this short 
time scale. Retamal and Zagury show that this re- 
mains valid under the Markovian and Born approxima- 
tions (Lindblad). At finite temperature a coherent state 
has a decoherence rate time scale given by 1/jt- We also 
see that, for this irreversible dissipation mechanism only 
the initial field coherence, (a^)(a), allows to increase the 
decoherence time scale. 



III. SPIN-BOSON MODEL 

Now let us consider a two-level system interacting 
with an external laser mode of frequency ujf and with 



a bosonic bath modeled by an infinite collection of quan- 
tized harmonic oscillators having frequencies The 
Hamiltonian which drives the unitary dynamics of the 
whole system, in the Hr = u>f(cr z + X^fc fyl&fc) rotating 
wave frame has the form [h = 1): 



H = Acr z +^2A k blb k + n(a + +a-) 

k 

+ J29k(b k a+ + bla^ , (7) 

k 



where A = Aq — ojf and Aq = u>i — Uo represent 
the energy difference between the upper state |1) and 
the lower state |0) of the two- level system. Operators 
a z , <j±, obey the standard SU(2) commutation relations. 
Afc = ujk — u>f, and gk's are the respective coupling con- 
stants for the dipolar interaction between the fcth mode 
and the two-level system, ft stands for the Rabi fre- 
quency which determines the coupling to the external 
classical field, bt and b k are the creation and the anni- 
hilation bosonic operators of the fcth mode, respectively. 
We assume the external field to be near resonance with 
the transition of the two-level system whereas the boson 
modes are considered to be far off resonance in such a 
way that A& 3> lor- 

Therefore, boson mediated transitions, described by 
the fourth term on the right hand side of J7J, can 
be strongly suppressed. Thus, the effective Hamilto- 
nian approximately describing the interaction process 
can be obtained from the Q Hamiltonian by using the 
method of Lie rotations | 2li |22|, namely applying to 
the Q Hamiltonian the unitary transformation: U = 

ex PEfc £ fc(frfc (J + _ b l a -)}> witn £ fe = ffj/Ak < 1. Ne- 
glecting terms of order higher than e k we obtain the fol- 
lowing H e f f effective Hamiltonian: 



H eff = Aa z + &kb\b k + Q (cr H 



<i 



k 



E 

k,k> 



9k9k- 



n 



bk + b) 



(8) 



where we have considered that Aq — u>k ~ Aq. Thus, 
off-resonant transition terms have been eliminated and 
in their stead appear three diagonal terms in the a z rep- 



resentation: a constant shift only for the upper level (a 
classical Stark shift); a random phase shift term simi- 
lar to that which is described in Ref. [l^; and a pure- 
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dephasing term (that depends on the intensity of the bath 
modes) which gives rise to a randomized phase between 
the two levels, but what is striking about it is that here 
the coupling constant includes the Rabi frequency of the 
field whereas, in the standard spin-boson modes hamilto- 
nian, an external field is not coupled at all. Thus, from 
first principles, an effective Hamiltonian for describing 
the dispersive interaction of a two-level system with a 
bosonic reservoir is obtained. 

The transformation of H e ff into the interaction pic- 
ture separates the motion generated by 



Hn = Ac 



(9) 



from the motion generated by the interaction term V = 
H e ff — Ho, so that 

V = e iHoty e -iH t = a Q + B(t)a z , (10) 

where a(t) = n(a+e lAt + <7-er lAt ), and 
B(t) = 2ttJ2 jr ( b i elUkt + bke^ kt ) 

k G 

+ J2^(bl bk ,e^-^-n k 5 k , k ,)(ll) 

k,k' 

with hk being the average boson number of the fcth mode. 
Here we emphasize that the effective interaction hamil- 
tonian 111) II has the same form of the ohmic spin-boson 
model |7lll4j. 

Considering Eqs. © and I|1U|) we obtain the following 
td time scale order for this model: 



1 



2J(B?) R (l-(^\a z \^) 



(12) 



where we have supposed that each boson mode is initially 
in a thermal state at absolute temperature T. Under that 
condition the operator (|11|) satisfies (B)r = 0, and 



(B 2 )r = ^^3fc(2n fe 

G k 



1) 



£#(' 11 »«■ 



k,k' 



A 2 

G 



For a fixed bath state the eigenstates of a z , that is, 
the poles of the Bloch sphere, are stable under decoher- 
ence effects. Meanwhile, the states on the equator of 
the Bloch sphere are the less stable. If there are no spin- 
boson interactions, that is g k = Vfc, all the \xjj) states are 
stable since the classical field drives a unitary evolution. 
The classical field affects decreasing the td time scale or- 
der through an effective dispersive interaction similar to 
the pure-dephasing mechanism. We can distinguish two 
regimes: the first one appears at weak field regime, this 
is, <C g k ; the second one turns up at strong field limit, 



51 ^> gk- At the high T temperature limit: the strong 
field regime has associated a given by: 



A 



G 



^ly/2kT^(l-{i/;\(T z \^) 
whereas in the weak classical field limit the td becomes 



td,€l<^g k 



2fcT 7 ^l - {^\<j z 



where 7 = J^ fe g?/u! k . Here we obtain an important dif- 
ference for the td behavior in the two regimes: in the 
strong field limit td is proportional to 1/y/T; meanwhile, 
in the weak field regime td is proportional to 1/T, both 
at high temperature limit. 



IV. CONCLUSIONS 

In summary we have found a general expression which 
allows us to estimate the decoherence time scale order for 
which the instability of a quantum state is onset. This 
decoherence time scale is a functional of the initial states 
of both interacting subsystems, and of the kind of in- 
teraction between them. It is worth pointing out that 
the found decoherence time scale is not symmetric with 
respect to the involved subsystems. 

Specifically, for a spin-boson weak and dispersive inter- 
action driven by a classical field we have found diverse 
behaviors of td with respect to the absolute temperature, 
depending on the intensity of the external classical field. 

One physical system which for instance fits these con- 
ditions would be a semiconductor quantum dot, where 
the uppermost valence band and the lowest conduction 
band can be represented by the ground and the excited 
eigenstates of a z . It is well known that the energy of a 
phonon is smaller than the transition energy of a two- 
level quantum dot system [23I 0|. Thus, the far from 
resonance constraint that we have imposed on the spin- 
boson model would be well satisfied in a semiconductor 
quantum dot. 

Further studies could involve other natural interactions 
and other atomic configurations as well as a resonant 
quantum driven field mode. Besides, one could study 
a more general model for the decoherence mechanism, 
which would include both dephasing and dissipation, that 
is, considering two infinite sets of modes, one set around 
the resonance and the other one far from the resonance. 
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